The LOPES experiment — recent results, status and perspectives 
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Abstract 



The LOPES experiment at the Karlsruhe Institute of Technology has been taking radio data in the frequency range from 40 to 80 
MHz in coincidence with the KASCADE-Grande air shower detector since 2003. Various experimental configurations have been 
— , employed to study aspects such as the energy scaling, geomagnetic dependence, lateral distribution, and polarization of the radio 
j> emission from cosmic rays. The high quality per-event air shower information provided by KASCADE-Grande has been the key to 
many of these studies and has even allowed us to perform detailed per-event comparisons with simulations of the radio emission. 
2jJ In this article, we give an overview of results obtained by LOPES, and present the status and perspectives of the ever-evolving 
experiment. 

0^ "Keywords: cosmic rays, extensive air showers, radio detection, interferometry 



1. Introduction 
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> ■ 

k>( Forty years after the initial discovery of radio emission from 
j_j cosmic ray air showers yj], cosmic ray radio detection has once 
again become a very active field of research. The LOPES ex- 
periment [2] in particular has revived the radio detection ac- 
tivities with an innovative approach combining digital radio- 
interferometry with detailed air shower measurements employ- 
ing a classical particle detector array. The experiment is situated 
at Campus North of the Karlsruhe Institute of Technology, at 
the site of the KASCADE-Grande air shower experiment yHH]. 
The close integration of LOPES with KASCADE-Grande, lead- 
ing in particular to the availability of high-quality per-event air 
shower parameters, has proven to be the key to many of the suc- 
cesses of LOPES. In comparison, the CODALEMA experiment 
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in Nancay is situated at an observing site with a much qui- 
eter radio background, facilitating the detection of radio signals 
in individual antennas without resorting to interferometric tech- 
niques. The CODALEMA experiment, however, does not have 
access to per-event air shower information of as high a quality 
as LOPES. 

In this article, we first describe the basic properties of the 
LOPES experiment, followed by an overview of the analysis 
procedures we apply to our data before we finally discuss the 
different phases and results of the ever-evolving experiment. 

2. Basic properties 

LOPES is a digital radio-interferometer measuring in the 
40-80 MHz frequency window. This frequency window was 
chosen to avoid FM radio-transmitters at frequencies above 
80 MHz and short-wave and atmospheric noise at frequencies 
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below 40 MHz. All LOPES channels are continuously sampled 
with 12 bit ADCs and a sampling rate of 80 MHz, i.e., in the 
second Nyquist zone. The data are stored in ring buffers, and 
when a trigger arrives from the KASCADE or Grande arrays, 
0.8 ms of data are read out and stored for each channel. Events 
with energies above w 10 16 eV are triggered; the threshold for 
radio detection lies significantly higher. While these basic prop- 
erties are universal to all configurations of LOPES, the layout 
and type of antennas has changed over the course of time. 

In Fig. Q~l an overview is given of the different phases of 
LOPES: LOPES 10, LOPES 30, LOPES 30 pol and LOPES 3D. 
The different phases were tailored to different scientific ques- 
tions, and will be discussed in some more detail in the following 
sections. 

Although the layout of the antenna array has been chang- 
ing over the different phases, in all cases the antennas have 
been placed in a region coinciding with the KASCADE array 
of KASCADE-Grande, as depicted for the LOPES 30 pol setup 
in Fig. [2] The antennas used for all phases except LOPES 3D 
have been inverted-V dipole antennas based, like most of the 
electronics used, on prototype designs for LOFAR |0]. Further 
technical details on the hardware of LOPES can be found in 
reference Q]. 
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Figure 2: Layout of the LOPES 30 pol setup. Triangles denote east-west- and 
north-south-aligned antennas, respectively; stars denote positions with both an 
east-west and north-south-aligned antenna. 



3. Analysis procedure 

The analysis procedure of LOPES data exhibits a high degree 
of sophistication, which is necessary in particular because of the 
comparatively high radio-frequency interference background at 
the "industrial" location of the KIT campus north. Here, we 
shortly discuss the analysis steps applied in our standard recon- 
struction. 



1. An absolute amplitude calibration is applied to the data. 
As of LOPES 30, this is based on a measurement with an 
external, calibrated reference source (8J]. The remaining 
systematic uncertainty of electric field amplitudes between 
events is of order 10 %. 

2. To suppress narrow-band transmitters, a digital filtering 
algorithm suppressing peaks in the frequency-spectra of 
each antenna is applied. This hardly affects the broad- 
band cosmic ray radio signals, and therefore increases the 
signal-to-noise ratio significantly. 

3. The delay of each individual channel (dominated by ca- 
ble delays and determined in a dedicated calibration cam- 
paign) is corrected for. 

4. A determination of the relative phases of the beacon trans- 
mitter operated in the vicinity of the LOPES antennas 
(cf. Fig. [2]) is used to constrain the relative timing of the 
LOPES channels to * 1 ns precision [9]. Such a high 
precision is a necessary prerequisite for a reliable inter- 
ferometric analysis. Before the availability of the beacon 
transmitter, a TV transmitter in the LOPES band was used 
for the phase calibration. 

5. Dispersion introduced by the instrumental response (in 
particular the bandpass filters) is de-convoluted from the 
signal traces, increasing the signal-to-noise ratio of cos- 
mic ray radio pulses. 

6. The 12.5 ns sampled raw data are up-sampled (i.e., inter- 
polated correctly) to a higher sampling rate. This is possi- 
ble because the data were sampled correctly in the second 
Nyquist zone and therefore contain the full information in 
the 40-80 MHz frequency window. 

7. A digital beam-forming is applied to arrange the time- 
series data of all channels correctly for radio emission 
received from the presumed arrival direction. (For the 
first iteration, the direction reconstructed by KASCADE 
or Grande is used as a starting point. Likewise, the core 
position is taken from the KASCADE or Grande recon- 
struction.) 

8. The antenna characteristics (i.e., frequency-dependent 
gain) for the established arrival direction is corrected for. 

9. The cross-correlation beam is calculated using eq. ([TJ. Af- 
terwards, it is block-averaged over 37.5 ns (approximately 
the width of the LOPES impulse response) and fitted with 
a Gaussian. The height of the Gaussian is used as a mea- 
sure for the amplitude of the overall coherent radio signal. 

10. For high signal-to-noise events, the maximum electric 
field amplitude in each antenna in a time-window around 
the CC-beam maximum is determined for the study of per- 
event lateral distributions. 

11. A correlation with a database listing events recorded dur- 
ing thunderstorms is performed. 

Steps 7-9 are performed in an iterative loop, during which 
the arrival direction of the emission and the spherical curvature 
of the electromagnetic front are varied to maximize the value of 
the CC-beam. Typical values for the spherical curvature of the 
electromagnetic wavefront are in the range from « 2-10 km J2]. 
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Figure 1: Time-line of the evolution of the LOPES experiment. LOPES STAR commenced in late 2005 and runs in parallel with the other phases of LOPES. 



A plane-wave assumption produces significantly poorer coher- 
ence. The CC-beam for each sample is calculated as 



cc[t] 



N, 



. N-l N 

— 2Z s, ' [ w ] 



p airs i=\ j>i 



(i) 



with N the number of antennas, s„[t] the time shifted (cf. step 
6 above) data of the denoted antenna, and TVpairs the number of 
unique pairs of antennas. The sign for each sample is the sign 
it had before taking the square root. 



4. LOPES 10 

The initial phase of LOPES consisted of 10 inverted- V dipole 
antennas oriented in the east-west direction. Only a relative 
amplitude calibration was available for these channels. With 
LOPES 10, the proof-of -principle for radio detection with a 
digital radio-interferometer was made [2]. Additional results 
were a near-linear correlation of the radio field strength with the 
primary particle energy, confirming the expected coherence of 
the radio emission, a clear correlation of the field strength with 
the angle between magnetic field and shower axis (the so-called 
"geomagnetic angle", both results published in |2[]), and an ap- 
proximately exponential damping of the radio emission with 
lateral distance with a slope parameter of Rq = (230 + 51) m 
[10]. In the latter publication, it was also demonstrated that ra- 
dio signals could be detected up to lateral distances of > 500 m 
and that as of 10 17 eV, LOPES has a high detection efficiency. 
Furthermore, it was demonstrated with LOPES 10 that very 
highly inclined air showers with zenith angles up to above 80° 
can be well detected with radio antennas 111 ill and that the radio 
emission from air showers can be significantly enhanced dur- 
ing thunderstorms, see Fig.|4]and [13]. In contrast, fair weather 
atmospheric electric fields do not significantly affect the radio 
emission. 



5. LOPES 30 

After the success of LOPES 10, the experiment was extended 
to 30 inverted-V dipole antennas aligned in the east-west direc- 
tion. The higher number of antennas distributed over a larger 
area increased the sensitivity (number of antennas) and angu- 
lar resolution (longest baseline) of the array. In addition, the 
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Figure 4: Comparison of events measured during thunderstorms (red circles) 
with a set of "twin" events with similar parameters taken during fair weather 
conditions (blue crosses). The y-axis quantifies the deviation of the measured 
field strength from the field strength expected for an event with the given pa- 
rameters. During thunderstorms, LOPES clearly sees enhanced radio emission 
from cosmic ray air showers. 



new layout was aimed at performing detailed per-event lateral- 
distribution studies. The reconfiguration of LOPES was com- 
plemented with regular calibration campaigns using an external 
reference source, thereby providing an unprecedented absolute 
calibration of each individual detection channel [8] with a re- 
maining systematic uncertainty between events of » 10% with 
respect to the electric field strength. To further improve the 
understanding of the instrumental response and possible envi- 
ronmental effects on the radio emission, a thorough monitoring 
of all relevant environmental parameters was started. This in- 
cluded in particular a monitoring of the atmospheric electric 
field to be able to reliably identify events recorded during thun- 
derstorms. 

Results from LOPES 30 are plentiful, and we only present a 
few highlights here. The correlations of the radio electric field 
strength with primary particle energy, lateral distance and geo- 
magnetic angle were studied once more using the absolute cali- 
brated, east-west polarization LOPES 30 data B12I1 . The results 
are depicted in Fig.[5]and can be parameterized with the follow- 
ing formula: 



st= (11 ± 1.) ((1.16 + 0.025) -cos a) cos 6 

cxpf -** )( E > f 95±004) r ^ ] 

U236+81)mA 10 I7 eV/ LmMHzJ 



(2) 
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Figure 3: Left: Up-sampled and beam-formed cosmic ray radio pulses seen in individual LOPES antennas. Right: Block-averaged cross-correlation-beam calculated 
from the individual antenna signals (blue) with its Gauss fit (light blue). The brown line denotes the power-beam, a quantity describing the total power received in 
all antennas together. As the cross-correlation-beam and power-beam coincide, the power in the radio signal is fully coherent. 



with a the geomagnetic angle, 6 the zenith angle, Rsa the mean 
distance of the antennas to the shower axis, and E p the primary 
particle energy. The given errors are the statistical errors from 
the fit. 

Another study dealt with the determination of the direction 
resolution, which corresponds to « 1.3° and is probably lim- 
ited by our current understanding of the shape of the electro- 
magnetic radio wave front 11 1 411 . Frequency spectra of cosmic 
ray radio signals were investigated with LOPES 30 as well 
which confirmed the expectation that the signal falls off 
to higher frequencies. 

A particular highlight of LOPES 30 results are the detailed, 
absolute calibrated per-event lateral distributions presented in 
fl^l . two of which are shown in Fig. [6] It was found that 
* 80% of events can be well-described with an exponential 
lateral distribution function with a typical exponential slope of 
Ro ~ 125 m. Up to » 20% of the events seem to show either 
completely flat lateral distributions or lateral distributions flat- 
tening towards the shower axis. The latter events occur predom- 
inantly at high zenith angles and for small lateral distances. In 
addition, it could be shown that a power-law parameterization is 
not able to describe the radio lateral distribution near the shower 
axis as well as an exponential distribution. These high-quality 
per-event lateral distributions are ideally suited for comparisons 
with radio emission simulations such as those performed with 
REAS2 Hzl] and REAS3 0, [H]. The agreement in particular 
with REAS3 simulations is remarkable with respect to both the 
absolute scale and the lateral slope, as illustrated by the exam- 
ple events shown in Fig.|7]and the histograms of Rq parameters 
for measurements and RE AS 3 simulations of both proton- and 
iron-induced air showers shown in Fig. [8] An exception are 
some of the measured events with very flat lateral distributions 
which cannot always be adequately reproduced by the simu- 
lations. Further systematic comparisons taking into account a 
more sophisticated treatment of detector effects are currently 
being performed. One effect that can influence such studies is 



the treatment of noise in the measured data, which has to be 
carried out with particular care Ir20ll . 

Another important question being studied with LOPES data 
is the search for a composition sensitivity of the radio emission. 
Such a sensitivity in various parameters accessible through ra- 
dio emission measurements has been predicted by simulation 
studies Ir2ll I22I1 . According to Ir2lll . the steepness of the ra- 
dio lateral distribution, in case of an exponential fit thus the Rq 
slope parameter, should be linked to the X m . dx of the air shower 
and hence the mass of the primary particle. KASCADE-Grande 
provides us with a reliable mass estimator given by the ratio of 
muon to electron numbers measured for each individual shower. 
The fact that air showers initiated by heavier nuclei exhibit a 
higher fraction of muons can be used in statistical analyses to 
decompose the air shower flux into different elemental groups 
Ir23ll : an identification on the basis of individual events is, how- 
ever, next-to-impossible due to shower-to-shower fluctuations. 
Preliminary investigations of a possible correlation of the radio 
lateral slope parameter Rq with the muon-to-electron-number 
ratio provided by KASCADE look promising (cf. Fig. [9]). A 
very detailed analysis of systematic effects will, however, be 
necessary before any definitive statements can be made. Fur- 
thermore, the prominence of a signature for radio emission sen- 
sitivity on the primary cosmic ray mass depends on the com- 
position of the cosmic rays in the energy range observable with 
LOPES (i.e., at >x 10 17 eV). If the cosmic ray flux at these 
energies is indeed iron-dominated, as can be expected from 
KASCADE-Grande results 112311 . the signature might be some- 
what weak. 



6. LOPES 30 pol 

After the measurements with 30 east-west-polarized anten- 
nas, LOPES was reconfigured to measure both the north-south 
and east-west polarization components of the radio emission 
from extensive air showers. The motivation for this reconfigura- 
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Figure 6: Examples of the radio lateral distribution of individual events measured with LOPES 30. The left example shows a typical event which can be reproduced 
well with an exponential function. The right example shows an event flattening towards the shower axis. Power-law parameterizations with a spectral index k 
describe the data close to the shower axis significantly worse than an exponential function. 




Figure 7: Comparison of two single-event LOPES 30 lateral distribution (blue) with REAS2 (green) and REAS3 (red) simulations of the same events. In the left 
plot, the REAS3 simulation is for a proton primary, in the right plot for an iron primary. The REAS simulations were based on the KASCADE reconstruction 
parameters and the simulations of "typical" showers with CORSIKA using QGSJET-II.03 and UrQMDl.3.1. 



tion was the fact that polarization studies are an ideal tool to dis- 
tinguish different radio emission mechanisms. LOPES 30 pol 
consisted of 15 channels with east-west polarization and 15 
channels with north-south polarization. At 5 antenna positions, 
both polarizations were combined at the same location. An 
overview of the LOPES 30 pol layout is shown in Fig. [2] Af- 
ter a TV transmitter in the LOPES frequency window which 
had been used as a reference to correct timing delays had been 
shut down, the beacon transmitter mentioned already in section 
[3] was set up to perform the phase correction with significantly 
improved accuracy |@]. 



First investigations of the data acquired with the 
LOPES 30 pol setup have confirmed that the polarization 
of the radio emission depends on the azimuth angle of the air 
shower arrival direction, see Fig. [10] and l24ll . Further studies 
are being carried out to analyze in detail the polarization 
characteristics of the radio emission and their relation to 
simplified v x § models or full-fledged simulations (see, e.g., 

lH). 



7. LOPES 3D 

In spring 2010, LOPES was once more reconfigured to the 
LOPES 3D setup. As the name implies, this setup is targeted at 
measuring the complete three-dimensional electric field vector 
instead of just a two-dimensional projection of the electric field 
on the horizontal plane. This has a number of benefits: 

• In contrast to a two-dimensional measurement, the com- 
plete radio signal is sampled, which leads to an increased 
signal-to-noise ratio for events with a significant fraction 
of vertical polarization, in particular highly inclined air 
showers. 

• As electromagnetic waves in the atmosphere are transverse 
waves, the determination of the complete electric field vec- 
tor provides information on the propagation direction of 
the emission already using a single antenna. This informa- 
tion can be used for cross-checks with the arrival direction 
determined by the arrival time distribution. 
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Figure 8: Histograms of the Rq distributions for a selection of LOPES 30 events in comparison with those of per-event REAS3 simulations for protons (left) and 
iron nuclei (right) as primary particles. The REAS3 simulations were performed on the basis of CORSIKA simulations using QGSJET-II.03 and UrQMD 1.3.1. 
The typical energy of the events is a 10 17 eV. 
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Figure 9: Correlation of the radio lateral slope parameter Rq with the ratio of muons and electrons measured by KASCADE for REAS3 simulations (left) and 
LOPES data (right). In the simulations plot, red points denote iron-initiated showers, blue points denote proton-initiated showers. The REAS3 simulations were 
performed on the basis of CORSIKA simulations using QGSIET-II.03 and UrQMD 1.3.1. The typical energy of the events is as lO 17 eV. 



• The complete information provided by the three- 
dimensional electric field vector can potentially help in the 
development of a polarization-based veto for an advanced 
radio self-trigger. 

• Emission models make predictions on the fraction of ver- 
tical polarization of the radio emission. In case of a sim- 
plified v x B model, e.g., the vertical component should 
be proportional to the north-south component. With 
three-dimensional measurements, such predictions can be 
tested. 

For the three-dimensional measurements, the antennas used in 
LOPES had to be replaced. The electronics, with the excep- 
tion of the low-noise amplifiers connected to the antennas, re- 
mained unchanged. Initial measurements with a SALLA an- 
tenna 12611 plus a vertical dipole and a self-developed tripole 
(see Fig. [iTIi consisting of three crossed dipoles have both been 
successful. For the final setup, the tripole antennas have been 



selected because they see a clearer signal of the reference bea- 
con and because they represent a more homogeneous design 
than a SALLA plus a vertical dipole, which helps to minimize 
systematic uncertainties related to the simulation of the antenna 
characteristics. Reconfiguration of LOPES to LOPES 3D began 
in February 2010. All calibration steps (determination of an- 
tenna positions, timing calibration, determination of reference 
phases, absolute amplitude calibration, simulation of antenna 
gain patterns, etc.) have been successfully performed. Since 
May 2010, the LOPES 3D setup with a layout of 10 tripoles as 
depicted in Fig.Q~2]is in stable data taking. 

8. LOPES STAR 

LOPES constitutes an ideal environment for research and de- 
velopment aimed at the large-scale application of the radio de- 
tection technique envisaged for example with the Auger En- 
gineering Radio Array (AERA) I127[ 12811 . Developments of 
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new antenna types for radio detection of cosmic rays such 
as the SALLA l26ll and a sophisticated digital self-trigger 
[29] have been carried out under the roof of LOPES with the 
LOPES STAR setup consisting of 10 dual-polarized antennas dis- 
tributed within the KASCADE-Grande array (see Fig. [2j>. 



9. Conclusions and Outlook 

The LOPES experiment has been very successful in study- 
ing the physics of radio emission from cosmic ray air showers 
with a variety of setups tailored to specific scientific questions. 
It has made key contributions to the understanding of the radio 




30 40 50 60 

Geomagnetic angle [degree] 

Figure 10: Measurements with the LOPES 30 pol setup show a distinct de- 
pendence of the relative strength of the north-south and east-west polarization 
components (here plotted as ratios of the corresponding CC-beam amplitudes) 
with both the azimuth angle of the air shower arrival direction (top) and the 
geomagnetic angle (bottom). 



emission physics and confirmed theoretical expectations such 
as the predominantly geomagnetic origin of the emission, the 
coherence of the radiation in the frequency band up to 80 MHz, 
the high angular resolution achievable with the radio technique 
and the reliability of radio measurements in all but the most 
extreme weather conditions. Furthermore, the absolute cali- 
bration of LOPES facilitates quantitative comparisons with the 
predictions of modern radio emission models, and indeed mod- 
ern models such as REAS3 are able to reproduce many of the 
features of the radio emission, including its absolute strength. 

Today, cosmic ray radio detection is in a phase of transition 
from first-generation, small-scale experiments to a full-fledged, 
stand-alone detection technique for cosmic ray air showers. 
Major challenges involved with this transition are a reliable and 
efficient self-triggering of radio signals, a precise determination 
of the primary energy from radio measurements, and the veri- 
fication of the predicted composition sensitivity. The future re- 
sults of LOPES will contribute to solving these challenges and 
continue to drive the development of the radio detection tech- 
nique towards its application on large scales for the study of 
high- and ultra-high-energy cosmic rays. 

Acknowledgements 

Part of this research has been supported by grant num- 
ber VH-NG-413 of the Helmholtz Association. LOPES and 
KASCADE-Grande have been supported by the German Fed- 
eral Ministry of Education and Research. KASCADE-Grande 
is partly supported by the MIUR and INAF of Italy, the Pol- 
ish Ministry of Science and Higher Education and by the Ro- 
manian Authority for Scientific Research CNCSIS-UEFISCSU 
(grant IDEI 1442/2008). 



7 




Figure 11: A tripole of the LOPES 3D setup. 
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Figure 12: Layout of the LOPES 3D setup consisting of 10 tripoles marked 
with blue stars. 
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